Quantum energy flow in atomic ions moving in magnetic fields 
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Using a combination of semiclassical and recently developed wave packet propagation techniques 
we find the quantum self-ionization process of highly excited ions moving in magnetic fields which 
has its origin in the energy transfer from the center of mass to the electronic motion. It obeys a 
time scale by orders of magnitude larger than the corresponding classical process. Importantly a 
quantum coherence phenomenon leading to the intermittent behaviour of the ionization signal is 
found and analyzed. Universal properties of the ionization process are established. 



Rydberg atoms in strong external fields represent an exciting and very active research area both experimentally as 
well as theoretically. During the eighties the main focus was on the hydrogen atom JO assuming an infinite nuclear 
mass which reduces the dimensionality of the system. However in general the atom possesses a nonvanishing center 
of mass (CM) motion in the magnetic field giving rise to a variety of two-body phenomena . Turning to charged 
two-body systems like, for example, the He + -\on the residual coupling of the CM and electronic motion is represented 
by an oscillating electric field term yielding five relevant degrees of freedom. One of the most striking effects caused 
by the two-body character of the ion is the recently found classical self-ionization process j| which occurs due to 
energy transfer from the CM to the electronic motion. Since it is well-known that quantization can severely change 
the effects observed in classical dynamics (see ref. and refs. therein) we develop in the present work a quantum 
approach to the moving ion in a magnetic field. 

We are interested in the regime of high level density, i.e. high-lying excitations, for both the collective as well as 
electronic motion which depend on a number of parameters (field strength, total energy consisiting of the initial CM 
and internal energies). The ab initio description of the quantum dynamics in the above regime goes even beyond 
modern computational possibilities and we thus seek a semiclassical approach that is capable of describing the essential 
physics of the problem. The total pseudomomentum ]5| is a conserved quantity associated with the CM motion. In 
spite of the fact that its components perpendicular to the magnetic field are not independent, i.e. do not commute, 
it can be used to find a suitable transformation of the Hamiltonian to a particularly simple and physically appealing 
form For the He + -ion it reads H — Hi + H2 + H3 where 
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Here m, Mo and M are the electron, nuclear and total mass, respectively, a = (Mo + 2m)/M and Q is the net charge of 
the ion. B is the magnetic field vector which is assumed to point along the z-axis. (R, P) and (r, p) are the canonical 
pairs for the CM and internal motion, respectively. The CM motion parallel to the magnetic field separates exactly 
and undergoes a free translational motion. 

Hi and H3 depend exclusively on the CM and electronic degrees of freedom, respectively. Hi describes the free 
motion of a CM pseudoparticle with charge Q and mass M. H3 stands for the electronic motion in the presence of 
paramagnetic, diamagnetic as well as Coulomb interactions which, in analogy to the hydrogen atom [Q, exhibits a 
variety of classical and quantum properties with changing parameters, i.e. energy and/or field strength. Hi contains 
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the coupling between the CM and electronic motion of the ion and represents a Stark term with a rapidly oscillating 
electric field 1/M (B x (P — Q/2B x R)) determined by the dynamics of the ion. This coupling term is responsible 
for the effects and phenomena discussed in the present investigation. 

The essential elements of our semiclassical approach are the following. Since we consider the case of a rapidly moving 
ion in a magnetic field a classical treatment of the CM motion coupled to the quantized electronic degrees of freedom 
seems appropriate: the CM is propagated with an effective Hamiltonian containing the corresponding expectation 
values with respect to the electronic quantum states. The latter obey a time-dependent Schrodinger equation which 
involves the classical CM trajectory. Both the electronic and CM motion have to be integrated simultaneously. The 
key idea of this semiclassical approach goes back to refs. where it has been applied to the dynamics of molecular 
processes. Our resulting time evolution equations read therefore as follows 

|p(*) = -^Wcj(R(*),P(*)) jR(t) = -^H cl (R(t),P(t)) a^(r,()=?< 5 (R(i),P(i),rWr,t) (4) 

with the effective Hamiltonian 

H c i(R, P) = Hi + (ip(r, t)\U 2 + n 3 \iP(r, t)) H q (R(t), P(t), r, p) = H 3 (r, p) + W 2 (R(t),P(t), r) (5) 

This scheme represents a balanced treatment of the coupled classical and quantum degrees of freedom of the ion and 
takes account of the energy flow among them. It possesses the important property of conserving the total energy 
which is particularly important for the correct description of the energy transfer processes occuring in our system. 
Since the typical energies associated with the fast heavy CM degrees of freedom are many orders of magnitude larger 
than the corresponding elementary quantum (TiQB/M) we expect the above scheme to yield reliable results. 

Our approach to the solution of the time-dependent Schrodinger equation, which yields the dynamics of an initially 
defined wave packet ip{r, t), is based on a recently developed nonperturbative hybrid method pO|-p^|. It uses a global 
basis on a subspace grid for the angular variables (9, (f) and a variable-step finite-difference approximation for the 
radial variable. The angular grid is obtained from the nodes of the corresponding Gaussian quadrature with respect 
to and </>, which is in the spirit of the discrete variable techniques yielding a diagonal representation for any local 
interaction fic[| . As a consequence one remains with the Schrddinger-type time-dependent radial equations coupled 
only through non-diagonal matrix elements of the kinetic energy operator. This vector equation is propagated using 
a splitting-up method [jl3|| , which permits a simple diagonalization procedure for the remaining non-diagonal part 
JTiyiSq . Our scheme is unconditionally stable, saves unitarity and has the same order of accuracy as the conventional 
Crank-Nickolson algorithm, i.e. ~ 0(At 2 ) where At is the time step size. In order to avoid reflections of the wave 
packet from the right edge of the radial grid we introduce absorbing boundary conditions. The extension of the radial 
grid is chosen (see below) such that it exceeds the center of the radial distribution of the initial wave packet by more 
than one order of magnitude. The typical frequencies associated with the motion of the Rydberg electron and the 
CM motion are different by several orders of magnitude (see ref. ||). To investigate the quantum energy transfer 
mechanisms requires therefore the integration of the above equations of motion for a typical time which is a multiple 
of the time scale of the heavy particle (CM) . This corresponds to many thousand cycles of the Rydberg electron. Such 
a detailed investigation would have been impossible without the use of the above-described combination of highly 
efficient techniques. 

We assume that the He + ion is accelerated up to some value Eqm of the kinetic energy of the CM motion and its 
electron is being excited to some Rydberg state nlm in field-free space. Thereafter it enters the magnetic field. In the 
following we choose Eqm = 100a. u., n — 25, 1 = rn = and a strong laboratory field B = 10~ 4 a.u.(23.5 Tesla). The 
intial CM velocity is vcm — 0. 1656a. u. and oriented along the a;— axis. We remark that taking the above values for 
nlm, B for the He + ion with the assumption of an infinitely heavy nucleus the corresponding classical phase space 
is dominated by chaotic trajectories. Figs, (la, lb) illustrate results for the propagation of the wave packet with 
increasing time. More precisely we show the intersection of the integrated quantity ty(p,z,t) = J \ip\ 2 d(j) along the 
cylindrical p-axis for z = (fig. la) and its intersection along the z-axis for p = (fig. lb). Fig. la demonstrates that 
the motion of the wave packet is confined by the diamagnetic interaction with respect to the p-direction, i.e. the 
direction perpendicular to the magnetic field. For any propagation time its value drops by several orders of magnitude 
at some outer value p c for the p coordinate. As we shall see below (see also figure 3) the variation of p c is accompanied 
by a corresponding change in the internal/CM energies thereby demonstrating the flow of energy between the CM 
and electronic degrees of freedom. Fig. lb demonstrates that there is for certain time intervals (see below) almost 
no decay of the wave packet for large distances of the z— coordinate. Therefore we encounter a significant flux of 
probability parallel to the external field. Reaching the boundary of our radial grid r m — 20000a. u. it is absorbed 
and considered to represent the ionized state. Having established the existence of an ionizing probability flux parallel 
to the magnetic field we immediately realize from fig. lb that this flux is by no means constant in time but varies 
strongly. To see this more explicitly and also to gain an idea of the overall decay of the wave packet we show in fig. 2 
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the decay of the norm of the wave packet for a time which roughly corresponds to one cyclotron period 2irM/(QB) of 
the free CM motion in the field. Fig. 2 shows apart from an overall monotonous decay of the norm, which is due to 
the quantum self-ionization process, an amazing new feature: the norm exhibits an alternating sequence of plateaus 
and phases of strong decay. The widths of the plateaus slightly increases with increasing time. This intermittent 
behaviour of the ionization signal from the moving ions is a pure quantum phenomenon, i.e. does not occur in the 
corresponding classical ionization rates || . Furthermore the calculation of the classical ionization rates for the same 
parameters (field strength, energies) yields a typical ionization time which is by two orders of magnitude smaller than 
the one obtained by the quantum calculation. The ionization process is therefore significantly slowed down through 
the quantization of the system which is in the spirit of the quantum localization processes shown to exist in a variety 
of different physical systems (see ref. Q and refs. therein). The observed slowing down of the ionization process 
represents one important difference of the classical and quantum behaviour of the moving ion which occurs in spite 
of the fact that we are dealing with a highly excited system. 

The obvious question of the origin of the intermittent occurrence of the plateaus and ionization bursts arises now. 
At this point it is helpful to consider the behaviour of the CM energy as a function of time which is illustrated in fig. 
3. Starting with Ecm = 100a. u. at t = we observe a fast drop of it for short times yielding a minimum of Ecm 
at approximately t\ = 4 x 10 7 a.u.. Thereafter it raises and reaches a maximum at approximately t 2 = 1.2 x 10 8 a.M. 
after which it drops again, i.e. it shows an overall oscillating behaviour. The "valleys" of Ecm coincide with the 
plateaus of the norm decay whereas the regions with higher CM velocities correspond to phases of a strong norm 
decay of the wave packet. The increase of the widths of the plateaus in the norm decay (see fig. 2) matches the 
corresponding decrease of the frequency of the oscillations of the CM energy. Since the total energy is conserved this 
clearly shows that the ionization bursts correspond to phases of relatively low internal energy (although certainly 
above the ionization threshold) whereas the phases of higher internal energy go along with the plateaus of the norm 
behaviour, i.e. the localization of the electronic motion. This provides the key for the understanding of the rich 
structure of the norm decay. The phase of high energy for the electronic motion means that the magnetic interaction 
strongly dominates over the Coulomb interaction. This makes the electronic motion approximately separable with 
respect to the motion perpendicular and parallel to the magnetic field. As a consequence the energy transfer process 
from the degrees of freedom perpendicular to those parallel to the external field are very weak i.e. the ionization 
process is strongly suppressed. This corresponds to an almost integrable situation for the ions dynamics. On the 
contrary for relatively low internal energies the Coulomb interaction is much more relevant and mediates together 
with the coupling Hamiltonian Ti. 2 the energy transfer from the CM to the electron motion parallel to the magnetic 
field. As a result we encounter a flow of probability in the +/-z-directions which corresponds to the ionization burst. 
During this period of motion a comparatively strong dephasing of the wave packet takes place. 

The intermittent behaviour of the ionization rate can therefore be seen as a quantum manifestation for the switching 
between different regimes of the internal energy corresponding to weaker or stronger Coulomb interaction. Pumping 
energy from the CM to the electronic motion weakens the Coulomb interaction and leads to the suppression of the 
ionization process whereas pushing the energy back to the CM motion decreases the internal energy and enhances 
the Coulomb interaction. To elucidate the time scale on which this process takes place we have computed the 
autocorrelation function C(t) =< ip(t)\ij)(Q) > where <> means integration over the electronic coordinates. As 
a result we observe a modulation and recurrence of the autocorrelation at a time scale t ~ 1.6 x 10 8 a.u. which 
corresponds approximately to the recurrence of the plateaus for the norm decay. The corresponding power spectrum 
shows a broad peak at a frequency uj « 3.5 x 10~ 8 a.u.. An important feature of the quantum self-ionization process is 
the approximate stability of the time intervals corresponding to the plateaus of the norm (no ionization signal) with 
respect to variations of the initial CM velocity of the ion. Our investigation shows that decreasing the CM energy 
from 100a. it. to 12a. u. leads to a decrease with respect to the distances between the plateaus, i.e. the difference of 
the norm values belonging to different plateaus, roughly by a factor of two. This corresponds to a significant slowing 
down of the ionization process. However the widths of the plateaus remain rather stable and represent therefore a 
universal quantity which is approximately independent of the CM velocity. Varying the field strength causes a change 
of both the distances between the plateaus and their widths. 

The quantum self-ionization process should have implications on the physics of atoms and plasmas occurring in 
a number of different circumstances. Apart from this it obviously suggests itself for a laboratory experiment (the 
lifetime of the Rydberg states exceeds the time scale of ionization by orders of magnitude) which should be very 
attractive due to the expected intermittent ionization signal which is a process revealing the intrinsic structure and 
dynamics of the system during its different phases of motion. 
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Figure Captions 

Figure 1: (a) The intersection ^(p, z = 0, t) = J \ip\ 2 d(f> along the p axis, (b) The intersection = 0, z, t) along 
the z axis. (Atomic units are used). 

Figure 2: The norm of the electronic wave packet as a function of time (in units of 10 8 atomic units). 

Figure 3: The CM energy as a function of time (in units of 10 8 atomic units). 
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